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A new approach to the parameterization of pion form factors is presented and for illustration applied to the pion 
vector form factor. It has the correct analytic structure, is at low energies consistent with recent high accuracy 
analyses of 7T7t scattering phase shifts and, at high energies, maps smoothly onto the well-known, successful isobar 
model. 
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£>. _ 1. Introduction 

o\ 

In recent years the knowledge about the low energy two-pion system has improved significantly, both 

experimentally as well as theoretically: for the low partial waves phase shift parameterizations of high 

accuracy exist from different dispersive analyses, either involving data only [1], or involving both data as well 

i— / ' as constraints from chiral symmetry [2] . The analyses are based on Roy or Roy- type equations that respect 

fvi . analyticity as well as crossing symmetry. Especially, left-hand cuts are included without approximation. 

In contrast to this, pion form factors or production reactions are often modeled either by sums of Breit- 
Wigners or improved versions thereof [3,4] or by the K-matrix formalism. In case of overlapping resonances 
unitarity gets violated by the former ansatz. The K-matrix provides a clear improvement compared to 
the Breit-Wigner parameterization, since two-body unitarity is built in. However, in general analyticity is 
violated. On the one hand, in the standard treatment not the full dispersive corrections are considered (in 
the expressions for the self energies only the imaginary parts and their analytic continuation are being kept 
and not the full expressions — c.f. Eq. (11) below), although some works include them (see, e.g., Ref. [5]), 
on the other hand, the left hand cuts are not treated properly — if they are included at all, in order to fit 
the scattering amplitudes, they are often in the same way included in the production amplitude although 
there the left hand cuts are different or, as in case of form factors, even absent. 

To be specific, in this work we focus on form factors and scattering with the goal to present simple 
formulas that allow for a data analysis that is consistent with analyticity and unitarity, however, without 
the necessity to solve dynamical equations. In addition, we present formulas that, by construction, in the 
low energy regime map smoothly and consistently onto what can be derived from the high accuracy analyses 
nrentioned above and thus for the scattering even include the proper left hand cuts. 

As an example and for demonstration we apply the formalism in this paper to the pion vector form factor, 
related to 7T7t scattering in the p-wave. The experimental situation for the p resonances beyond the p(770) 
is at present not very clear: different experiments find indications for different resonances — for a summary 
of the current situation see 'note on the p(1450) and the p(1700)' in the Review of Particle Physics [6]. 
The formalism presented here could be an important step forward to clarify the situation for it allows for 
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Fig. 1. Diagrammatic representation for the various ingredients of the formalism. 

a simultaneous, consistent analysis of various channels/observables. To parameterize the vector form factor 
beyond s — 1 GeV 2 all studies agree on the need to include at least two resonances in addition to the /o(770), 
which is elastic. Thus, if we include one inelastic channel the pion vector form factor is parameterized in 
terms of in total 9 parameters — 8 for the resonances and 1 additional parameter for the p-u mixing (see 
Sec. 3). The number of parameters needed for each resonance agrees to standard parameterizations. The 
advantage of the parameterization presented here is that we do not need to approximate the left hand cuts, 
the consistency with low energy phase shift is ensured by construction and the connection between scattering 
and form factors is properly implemented. We found that once simultaneously to the form factor also the 
data on e + e~ -^non 2ir in the isovector state was fitted, the inclusion of not only an additional resonance 
but also of the ttuj cahnnel together with a direct tttt — » ttlo coupling became necessary. This model, with in 
total 17 parameters, allowed for an acceptable fit to the data. 

The paper is structured as follows: the most important formulas are motivated and presented in Sec. 2. 
Their more detailed derivation is given in Sec. 3. Results are presented in Sec. 4 and the paper closes with a 
short summary in Sec. 5. The two-potential formalism which forms the basis for the derivation is introduced 
in the Appendix. 

2. Summary of most important results 

When it comes to the interplay of resonance contributions and background terms the case of only one 
scattering channel is especially simple for here it is possible to give a closed form expression for the form 
factor T solely in terms of the elastic scattering phase shift — the so-called Omnes solution [7] . It is derived 
from a dispersion relation using the fact that 

disc(J"(s)) = 2iaT(s)* J"(s) , (1) 



where T denotes the on-shell elastic scattering amplitude, and a = y/1 — Am 2 / s the two-body phase space — 
for simplicity we assume the scattering particles to have equal mass m. Here 'disc' denotes the discontinuity 
of the form factor defined via 

disc(J"(s)) = T(s + ie) - T{s - ie) = 2Tm(J r ) . 

If the elastic phase shift is 5(s), with 

T(s) = -sm(6(s))e i5{s) , (2) 

a 

then the form factor reads in the absence of bound states 

T{s)=i{s)^\{ S )P A { S ) (3) 

with the Omnes function 

W .Ji/"*M}, (4, 

I 7T J 4m 2 S' S'-S) 

where the centrifugal barrier is introduced via the factor £(s) = y/s — s t hr for L-waves, where Sthr denotes 
the location of the threshold and the function Pa(s) is a polynomial. Its degree may be fixed by the large s 
behavior of the form factors. For values of s < 1 GeV 2 this methodology was used by various authors for the 
pion vector form factor, see, e.g. Refs. [4,8-13]. Here we will present a formalism that provides expressions 
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Fig. 2. Fits result for the pion p— wave phase shift (left panel) and inelasticity (right panel). The red solid (green dashed) line 
denotes the result of the fit #2 (#1). The dot-dashed line in the left panel refers to the input phase 5. Data are from Ref. [14] 

(solid dots — only data below 1.4 GeV are shown [15]), Ref. [16] (solid squares for solution ( ); solid triangles for solution 

( 1 )), and Ref. [17] (open dots) and Ref. [1] (turquoise band). 

that smoothly map onto Eq. (4) at lower energies while being an analytically improved version of the isobar 
model at higher energies. 

Eqs. (1) and (4) apply only if the interactions are purely elastic — for the latter it is even necessary that 
they are elastic up to infinite energies. Clearly this is not realistic. However, experimental data show that 
at higher energies inelasticities are typically accompanied by resonances. We therefore split the full, partial 
wave projected, interaction potential V into two pieces 



V{s) l0 = V{a)ii + V R (s) 



ij i 



(5) 



where i and j denote the channels. The crucial feature for this approach is that the potential V needs to 
be specified at no point. All what is needed are the corresponding phase shifts S. We now postulate the 
following properties: 

• the potential V is purely elastic, such that V is non-vanishing only for i = j = 1; 

• deviations in the tttt phase-shifts from S come either from s-channel resonances or via the coupling to 
inelastic channels; 

• all long ranged forces are in the elastic tttt interactions V; all interactions in other channels are regarded 
as short ranged. 

These are the model dependent assumptions of this approach (note, in case of the traditional isobar model 
one needs to assume that all interactions are mediated by s-channel resonances). Based on these Eqs. (1) 
and (4) can be easily generalized to multiple channels providing a convenient parameterization for both 
scattering as well as production amplitudes. 

The full scattering T-matrix appears as the solution of a Bethe-Salpeter equation with input potential 
V defined in Eq. (5). Using the two potential formalism (see Appendix) it is straightforward to derive the 
decomposition 



T{s)ij = 5ij5uT{s) + T R (s)ij = Sij5iiT(s) + r out (s)itij(s)yT in (s) 



(6) 



where T is the purely elastic scattering T-matrix that derives from the potential V. For the vertex 
functions one has from time reversal invariance r out = T[ n . They are, by assumption, non-trivial only in the 
tttt channel, where elastic scattering is mediated by T. As a result r out is diagonal. We may therefore write 
for any given partial wave 



f = —e fs sm(S) (7) 

and 

r ou t(s)» = &( S M<$]( S ) f° r i= l (^-channel); otherwise r out (s)j = &(s) , (8) 

where £i(s)y/s — s thr« is the multichannel version of the centrifugal barrier factor given above. The reso- 
nance T matrix tn may be written as 

*«(«)« = [1 - Vnis^s)]-,} V R (s) kj (9) 

with the resonance potential (note: not all resonances are in Vr; elastic resonances may be included in T 
as the p(770) in the example below) 

n (l) (l) 

Vb(*)« = - S 9% _ 9 \ ! Vn(«) - Vb(«) - Vb(0) + 7j*i + 7i*ji • (io) 

j=i s (0 

In Fig. 1 a graphic representation of the various quantities is given. The potential is subtracted at s = to 
ensure that the phase of the full T matrix at low energies agrees to the input phase S. Clearly, the procedure 
does not guarantee a priori that the phase of the full T-matrix is close to that of T in the whole range where 
T is well determined, however, in practice this is indeed the case: as can be seen from the left panel of Fig. 2, 
at energies below 1 GeV all curves shown, including that for the input phase shift, are indistinguishable. 

The term containing the parameter 7j, with "fi — 0, is included to allow for a direct transition of the irir 
channel (channel number 1) to inelastic channels. For the pion vector form factor the fit to data requires 
such a coupling to the ttlu channel in effect it provides a direct piruj coupling, which from phenomenology 
is known to be significant [18]. 

The self energy loop Ej in channel i appearing in Eq. (9) and Fig. 1 may be expressed via a twice subtracted 
dispersion relation (see Sec. 3 and Appendix) — one finds 



s i(s) = — — 



2 ^ds'a^OIW)! 2 



le 



(11) 



For given interactions Ty it is straight forward to calculate the form factors T{. As shown in the next 
section, one finds 

T{$)i - T out ( s ), [1 - Vnis^is)}^ 1 M k , (12) 

where the Mi denote point like source terms for the production of particles into channel i. It may be 
written as 

"-<*- Erbs-- < 13) 

i=i (0 

Since we assume that all interactions but those in the two-pion channel are driven by resonances, we choose 
cj. = cSik, with some constant c. In case of the pion vector form factor, discussed in detail below, charge 
conservation demands c = 1. In addition, to ensure the proper normalization of the form factor and to 
suppress the influence of the higher resonances on low energies, we use a photon-resonance coupling linear 
in s 1 . In practice this means replacing ay> in Eq. (13) by set"'. The ingredients of the formalism are 
illustrated in Fig. 1. 

All parameters introduced — k, g\ , a"', rnm, 7j — are real as long as all (relevant) channels are treated 
explicitly as a consequence of time reversal invariance. Thus, in case of two channels, for each inelastic 
resonance we need to include 4 real parameters. Each additional channel adds in one more parameter per 
resonance. For the predominantly elastic resonances, which are included in the vertex functions T, the only 
free parameters are their couplings to the inelastic channels introduced via the 7$. 

Eq. (12) is the central result of our paper. If all vertex functions were chosen to be constant, it would 
reduce to the famous P-vector formalism [19]. However, since in our case Ti, which enters explicitly in the 



1 On the Lagrangian level this means a coupling of the resonance to the photon field via F^d^Vv with F^ v for the electro- 
magnetic field strength tensor and V v for the resonance field. 
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Fig. 3. Real and imaginary parts for the self energies Ej(s) defined in Eq. (11). Left panel: self energy for the elastic channel, 
Ei(s) — the solid (dashed) line shows the imaginary (real) part. Right panel: self energies for the inelastic channels. The solid 
(dashed) line shows the imaginary (real) part for the An channel, while the dashed— dotted (double dashed-dotted) line shows 
the imaginary (real) part for the um channel. 

expression for J-\ as well as through Ei, is non-trivial, Eq. (12) provides a generalization to the conventional 
treatment. 

3. Derivation of the Formalism 



In this section the expressions presented in the previous section are derived. It is based on the two- 
potential formalism used, e.g., in Ref. [20] to control the Coulomb- nuclear interference and rederived in 
the notation used here in the Appendix. The important aspect to observe is that, although the derivation 
starts from a potential V, introduced in Eq. (5), the final expressions can all be expressed in terms of the 
corresponding scattering phase 5. 

The two-potential formalism splits the T-matrix into two pieces, the elastic T— matrix T and the resonance 
T-matrix Tr (c.f. Eq. (6)). The latter quantity is the solution of a Bethe-Salpeter equation derived from 
the resonance potential Vr, however, here the intermediate state free propagators, used in the standard 
treatment, are to be replaced by interacting propagators. This is achieved via inclusion of the vertex functions 
r on the external legs as well as a proper modification of the intermediate meson loops connecting two 
insertions of the resonance potential — the so-called resonance self-energies Ej. It is straightforward to 
show that the imaginary parts of the self-energies read (see Appendix) 

Im(E i (*))=(T < |r i (*)| a . (14) 

Thus the self energies can be calculated from a properly subtracted dispersion integral. Since E(so) and 
<9E/<9s(so) at some so can be absorbed into the resonance masses and wave function renormalization con- 
stants (or, more appropriate for the present context: in the effective coupling constants of the resonances) 
we here use a twice subtracted version — c.f. Eq. (11). The self energy in the 2-ir channel needs as input 
only the vertex function Ti(s) discussed in the previous paragraph. The resulting Ei(s) is shown in the left 
panel Fig. 3. 

To parameterize the inelastic channels we use for i = 2 a structureless 47r channel via the phase space factor 
(72 = y/l — 16m 2 /s , which provides the proper scaling of the four-body phase space near the threshold. For 
the barrier factor, appearing in Eqs. (8), we use £2 = ys — 16m 2 . The resulting self energy ^(s) is shown 
in the right panel Fig. 3. Although the KK channel contributes significantly to isoscalar rnr interactions, 
it gives negligible contributions in the isovector channel [21]. Therefore in fit #2 we include as additional 
inelastic channel (i — 3) the irui channel. For a discussion on the possible role of the ttoj channel on the ttit 
inelasticity see Ref. [22,23]. We take 03 = A(s, m^m^) 1 / 2 / 5 an d £3 = A( s , m2 J ,™ 2 -) 1 ^ 2 - The resulting self 
energy Ea(s) is also shown in the right panel Fig. 3. The linear rise of the imaginary part of both E2 and 
E 3 comes from the centrifugal barrier terms £2 and £3, respectively. 



To calculate the form factor we may write (using the notation of the Appendix) 

T = M + TGM , (15) 

where Gi is the operator representation for the integration over all intermediate n-body states of channel i 
and the production vertices Mj were defined in Eq. (13). Inserting Eq. (6) into Eq. (15), we get 

T = M + TGM + T R GM = r out (l + t R T\ n G)M . 

To proceed we may use the definition of the self energy Sj = ± ini Gi, to write 

t K S = [1 - VrZ]- 1 VrT. = -1 + [1 - VrZ]- 1 . 

Here we needed to assume that the range of interactions in the production vertex and in the vertex functions 
of the resonances is similar in all channels, for only then the same loop integral Sj can be used as self energy 
contribution for the resonances as well as convolution integral of Mi and the resonance potential. From this 
we get 

T = r out [i-v fl s]- 1 M , 

which agrees to Eq. (12). 

It is important to observe that the expression given in Eq. (12) is consistent with the coupled channel 
version of the unitarity relation for form factors, Eq. (1), since 

disc(j-) = disc(r out ) [i - VrZ]- 1 m + r* ut disc([i - VrSJ-^m 

= 2if vr out [i - VrY]- 1 m + r: ut [i - Vr^t 1 Vfldisc(s) [1 - VrT,}' 1 m 

= 2i (f* + r* out [i - VrZT 1 v R r* out ) <jr out [i - VrZ}- 1 m (16) 



where in the intermediate step the unitarity relation for the vertex function, Eq. (1), and the self energy, 
Eq. (14), were used. 

An interesting observable is the ratio r of the total cross section for e + e~ annihilation into hadronic states 
with 1 = 1 other than 7r + 7r~ over cr e + e -_ i . 7r + 7r - — a compilation of this quantity can be found in Rcf. [21]. 
In this ratio the unitarization effects in the resonance T-matrix cancel largely. For example, in case of only 
one inelastic channel we get 

(17) 



/CT 2 r out2 \ {i-v r11 ^ 1 )m 2 + v r12 ^ 1 m 1 



Wr out i/ (i - Vr2 2 s 2 )Mi + Vr 12 z 2 m 2 

clearly being very directly sensitive to the resonance parameters. 

The vector form factor in the two pion channel is directly accessible from two reactions, namely from 
e + e-annihilation and from r decays. In the former case, in addition to what was discussed so far, the isospin 
violating mechanism of p-uj mixing needs to be included. It is visible as a striking narrow structure close to 
the peak of the corresponding form factor (c.f. inlay in Fig. 4). The inclusion of this mixing in the present 
formalism is straight forward — we here use a slightly modified version to what is used in Ref. [11], namely, 
for the neutral, ir + ir~ , channel 

F 1 (s)—^F 1 ( 8 )(l + K 2 S ). (18) 

We here use m u = 0.7826 GcV and T u = 0.0085 GeV. The strength parameter k is part of the fit. The 
only difference between the form factors in the neutral (from e + e~ annihilation) and in the charged (from 
r decays) channel is that in the latter we choose k = 0. 

4. Results for the pion phases, inelasticities and form factors in the p-wave 

To be specific we now focus on the pion vector form factor. The resonance parameters will be determined 
by a fit to data on the pion vector form factor as well as inclusive data on inelastic channels. In addition we 
need as input the elastic phase shifts 8, which largely fix the properties of the p(770). 
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Fig. 4. Fit result for the pion vector form factor Fy = T\/£\. Left panel: For the neutral channel. The red solid (green dashed) 
line denotes the result of fit #1 (#2). The black dot-dashed line shows the form factor derived from the Omnes function only. 
Data are from the reaction e + e~~ — > 7T+7T - presented in Refs. [24-26]. Right panel: For the charged channel. Data are from 
Belle [27] and CLEO [28]. 
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Table 1 

Parameter values for the mixing parameter k and the resonance parameters for the two fits performed. The uncertainties listed 
refer to the statistical uncertainty of the fit only. Masses and 73 are given in GeV and GeV -2 , respectively; all other couplings 
are dimensionless. 

For our work we will use for energies below s cut = 1.4 2 GeV 2 the central values for the phase shift provided 
in Ref. [1] - - see Eqs. (A7) and (A8) therein. For energies above this value we smoothly extrapolate the 
phase shift to a value of 7r via 

5(a) = tt + fc ut ) - it) ( ^7 ) ■ (19) 

It turns out that for A > 2 GeV the results are basically insensitive to the actual value used 2 . We thus 
chose A = 10 GeV in what follows. The asymptotic value 7r for the phase ensures that the vertex function 
Ti(s) decreases as 1/s as demanded for the vector form factor. The resulting elastic tttt phase shifts shown 
as the black dot-dashed line in the left panel of Fig. 2. The form factor from the Omnes function alone is 
shown as the black dot-dashed line in Figs. 4. It provides an acceptable description of the data up to s = 1 
GeV 2 , although there are some deviations visible (see inlay in the left panel) . At higher energies a significant 
deviation becomes visible. 

Two different fits were performed using the MINUIT package of the CERN library: the goal of fit #1 is 
to find the minimal parameter set necessary to get an acceptable description of the pion vector form factor. 
In line with previous studies we find that at least two resonances and one inelastic channel are needed. In 
addition we may chose all % = (c.f. Eq. (10)). The fit result is shown as the dashed green line in Fig. 4 



2 This is correct up to an unphysical pole located at s 
s > -A 2 . 
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Fig. 5. Results for the ratio r = a 
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0.6 GeV 2 results from p-u mixing present in the 



for the pion vector form factor and in Fig. 2 for p-wave phases and inelasticity. However, it turned out that 
with this parameter set it was impossible to simultaneously also describe r — see green dashed line in Fig. 5. 
Therefore we performed a second fit, fit #2, with the goal to get a decent description of both Fy as well 
as r. This called for an inclusion of 2 inelastic channels — in addition to the 4.7T channel already included 
in fit #1 we now also include the 7rw channel — 3 resonances as well as a non- vanishing value of 73. Fit 
#2 has 17 adjustable parameters. The fit results are shown as the red solid lines in Figs. 2, 4 and 5 for 
phases and inelasticities, the pion vector form factor from e + e~ annihilation and r decays, and the ratio r, 
respectively. Note, in case of the pion vector form factor we fit to the BaBar data on e + e~ annihilation only, 
for it extends to higher energies 3 . Especially, the r data are not included in the fit — the result comes out 
as a prediction. The parameters determined in both fits are given in Tab. 1. It is important to note that for 
the pion vector form factor alone both fits are of similar quality: for fit $1 and #2 we have \ /d.o.f=1.3 and 
1.2, respectively. However, only fit #2 provides an acceptable description for r. This result nicely illustrates 
that one should not analyze the pion vector form factor without looking at the non 2ir channels at the same 
time — to our knowledge in this work a combined analysis was performed for the first time. Note that at 
s = 2 GeV 2 Ref. [21] reports a value of r ~ 26, which shows that already at this relatively low energy the 
2tt channel provides only a small fraction of the e + e _ annihilation rate in the I — \ channel. 

The s dependence of the non 27r data shown in Fig. 5 calls at least for two inelastic channels, since there is 
a change in slope visible at around s ~ 0.9. In fit #2 this is accounted for by the inclusion of the ttuj channel. 
However, even this three channel fit is still too simplified, for there should be not only correlations amongst 
the 4 pions in channel 2 included, e.g. from ai7r, pp, and pa. In addition there are also channels like r]mr. 
The data included in the current study does not allow one to disentangle these and therefore to improve the 
description of r further. What is necessary is an inclusion of the large number of exclusive measurements 
available from e + e~ annihilation. We leave this to a future study. 



3 In this exploratory study we regard this as appropriate since in this work we do not perform an uncertainty estimate of the 
parameters extracted. 
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Fig. 6. Left panel: Fit result for the pion vector form factor at space-like energies. The lines are the same as in Fig. 4. Data 
are from Ref. [29,26]. Right panel: phase difference between ip, the phase of the form factor, and the scattering phase shift 8. 
The data shown indicate the upper bound of the phase shift difference presented in Ref. [21]. 

It is important to stress that the mass parameters given in the table are bare parameters that get renor- 
malized by the self energies. It is therefore possible that the resonance poles come out similar in both fits, 
for the linearization effects are different. However, we postpone the determination of pole positions and 
residues, which requires the evaluation of the elastic T-matrix T in the complex plane, to a later work. 

At small, space like energies the results of the two fits give identical results, c.f. left panel of Fig. 6. Both 
lead to a pion radius slightly enhanced compared to what comes from the Omnes function itself, shown by 
the dot-dashed line: through the inclusion of the high-lying resonances the mean square charge radius of 
the pion increases by nearly 10% from 0.40 fm 2 to 0.44 fm 2 . The latter value is consistent with the values 
extracted in Refs. [30,31]. For the curvature we find cy = 3.9 GeV -4 , in line with Refs. [13,32]. The effect of 
the higher resonances is quantitatively in line with expectations from dimensional analysis that predicts an 
effect on the mean square radius of order of the square of the inverse resonance mass ~ 0.02 fm 2 . At higher 
space like energies the results move appart from each other. While fit # 2 is consistent with the largely 
model independent bounds for the form factor derived in Ref. [33] , fit # 1 (the form factor from the Omnes 
function) is at most marginally consistent (inconsistent). 

The left panel of Fig. 2 nicely illustrates that within the formalism presented the high accuracy phase 
shifts up to 1 GeV are reproduced very well. One also finds that the phase shifts for the full model largely 
agree to the input phase in the whole energy range considered (this is not the case for the phase of the form 
factor, as discussed in the next paragraph) as well as to the data of Ref. [14,16] 4 . This happens, since 
the resonance couplings to the tttt channel are rather small — the resonances show up prominently in the 
form factor only due to large couplings to the photon. Especially, the present model can not account for the 
significant inelasticity visible in the data of Ref. [16] — the according to Ref. [34] preferred solution ( — | — ) 
shows an inelasticity of the 0.8 already at 1 GeV. At this point in time it is not possible to decide whether 
this failure is an indication of a short coming of the model used here, or of the data of Refs. [14,16]. What 
might support the latter conjecture is that the data on r\ of Refs. [14,16] are in disagreement with both the 
analysis of Ref. [17] as well as that of Ref. [1] at around 1 GeV. 

In the elastic regime the phase of the form factor has to agree to the phase of elastic scattering — a fact 
known as Watson theorem. At higher energies this connection is lost. In the right panel of Fig. 6 we show the 
difference between -0, the phase of the form factor, and the scattering phase shift S. Also shown in the panel 
is the allowed upper bound of the phase shift difference given in Ref. [21]. As one can see our amplitudes 



Only two of the 4 solutions presented in that paper are given, 
shifts of Ref. [1] . 



for the other two are in strong discrepancy with the phase 



largely exhaust the range allowed by unitarity. This reflects again the fact that in the present formalism all 
resonances besides the p(770) couple to elastic scattering only weakly. 

5. Summary and Outlook 

In this paper a formalism was presented that allows for a simultaneous description of both irir scattering 
data as well as form factors without the need to model the low energy regime: at low energies tttt phases can 
be used as input directly. At higher energies the formalism maps smoothly onto the well known N/D method 
which is similar to the iiT-matrix approach, however, with improved analytical properties. As an example in 
this paper the formalism was applied to pion pairs in the p-wave. An excellent description is found for the 
pion vector form factor in both the neutral channel — from e + e~ annihilations, with p-uj mixing included 
— as well as the charged channel — from r decays. In addition we also found a qualitative agreement with 
data on the non 2n channels from e + e~ annihilations — these data were studied within a dynamical model 
here for the first time. 

We found, however, that within the given formalism it was not possible to to describe the behavior of the 
inelasticity given in Refs. [14,16]. At this point in time we are not able to judge if this deviation indicates 
a short coming of the model or points at a problem in the data. However, the observation that the values 
of rj of Refs. [14,16] arc in disagreement with the analyses of Refs. [1,17] at s ~ 1 GeV 2 might indicate that 
there is a problem in the data of Refs. [14,16] also at higher energies. 

The formalism described here can be applied to all partial waves, especially also the isoscalar s-wave. 
Here, however, it is less clear what to use for the elastic phase shift S, since the pronounced structure from 
the /o(980), which also couples strongly to KK, shows up already short after the phase reached 90°. We 
leave this study to a future work. 
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Appendix A. The two potential formalism 

Let us assume that there is a sensible way to split the scattering potential into two pieces (as in the main 
text the potentials V and V R , the vertex function T as well as the T-matrices are matrices in channel space, 
while the form factor F and the production vertices M are vectors in channel space) 

V = V + V R . 

We will show in this Appendix that the full T-matrix can be split accordingly. In operator form the Bethc 
Salpeter equation for the T matrix may be written as 

T = V + VGT = V + V R + (V + V R )GT . 

Here Gi denotes the operator for the integral over the n-particle intermediate state of channel i, e.g. for the 
two-7r intermediate state we have 

1 f d A k 1 1 

VGVcc -J J^ V ^ ■■ ) ki- m 2 + ie{k-Py-mi + ie V{K - ] ' 

where P denotes the total 4-momcntum of the system, P 2 = s. Note, not all arguments of the potential V 
are shown explicitly. Introducing T as the solution of 

f = V + VGT 

and the dressed vertex functions (to simplify notations, in the appendix we do not show the centrifugal 
barrier factors £ explicitly) 

Tout = 1 + TG and rf n = 1 + GT 
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we get 

T = f + T R = f + V R Tl + (V R G + VG)T R . 

Due to time reversal invariance we have r out = T\ n . Since disc(G)=2ia and disc(T) = 2iaT*T one has 

disc(r out (s)) = disc(f )G + f *disc(G) = af(s)*T out {s) , 

thus r out holds the unitarity relation for a form factor in a channel where the interactions are given by T . 
We may therefore dehne 

Tr = r out iflr in 

and derive with E = GT out and VGT out — T out — 1 

t R = v R + v R n R — ► t R = [i - VrT,}- 1 v r . 

From the definition above the discontinuity of the self-energy £ is found to be 

disc(S) = disc(G)r out + G*disc(r out ) = 2t (l + G*f *) aT out , 



which was used to derive Eq. (11). 
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